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Abstract. The strange quark matter (SQM) may be the ground state of 
QCD; nuggets of SQM could be present in cosmic rays (CR). SLIM is a large area 
experiment, using CR39 and Makrofol track etch detectors, presently deployed 
' at the high altitude CR Laboratory of Chacaltaya, Bolivia. We discuss the 

, expected properties of SQM, from the point of view of its search with SLIM. 

■» ' We present also some preliminary results from SLIM. 
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SLIM is a large area experiment (440 m^) installed at the Chacaltaya CR lab 
since 2001; an additional 100 ir? were installed at Koksil, Pakistan, since 2003 ^ 
1, 31 El El- With an average exposure time of about 4 years, SLIM would be 
' sensitive to a flux of downgoing exotic particles at a level of 10"-'^^ cm~^sr~^s~^. 

The main goal of SLIM is the search for intermediate and low mass magnetic 
monopoles, but it can be sensible to other exotica, as nuggets of SQM (known 
as "nuclearites" or "strangelets" ) , Q-balls, etc. 

We focus on the search for nuggets of SQM (both in the low mass region, 
in which they are expected to behave more or less like super-heavy nuclei, as 
well as in the intermediate mass region, that is if they are heavy enough to 
behave like micro-meteorites, but not so large to be able to cross the Earth) 
in the cosmic radiation, using the SLIM detector. After a short description of 
the experiment and of the analysis procedures, we briefly review the expected 
properties of SQM, relevant from the point of view of SLIM. We then investigate 
the arrival and detection conditions of cosmic ray SQM nuggets in SLIM and 
discuss some preliminary results yielded by the analysis of a part of the detector. 

^This paper refers to the Chacaltaya location only. 
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Figure 1: A sketch of a SLIM module. The stars indicate the detectors in which 
an etchable track is prodviced by the passage of different particles. 

2 The experiment 

SLIM (Search for Light and Intermediate mass Monopoles) is a large area de- 
tector (about 440 m^) presently exposed to penetrating CR at the Laboratory 
of Chacaltaya (Bolivia), at an altitude of 5230 m a.s.l. The detector consists on 
modules of 24 x 24 cm^, made of 3 sheets of CR39 and Makrofol nuclear track 
detectors (NTDs), and an aluminium absorber. This stack structure is similar 
to that used by the MACRO track-etch subdetector [51 [7j . Fig. 1 shows the 
sketch of one module. Two additional Lexan sheets are placed on the top and 
the bottom of each stack; those detectors are not used in our analysis, their role 
is only to absorb some of the a particles produced by the Radon decays in the air 
around SLIM. The response of NTDs depends on the environmental conditions 
during their exposure; thus each module is sealed in a mylar bag, filled with dry 
air at the normal atmospheric pressure (note that the atmospheric pressure at 
Chacaltaya is 0.5 atm.). 

The working principle of NTDs is qualitatively presented in Fig. 2. When a 
particle crosses such a detector, if the restricted energy loss (REL) is above some 
specific threshold, the polymeric structure is affected along is path yielding the 
so called "latent track" . In the particular case of a particle of electric charge Z 
and velocity (5 — vjc, the produced damage is a function of Zj i3. The tracks 
became visible after chemical etching, typically in aqueous solutions of NaOH 
or KOH, as the etching velocity along the latent track (wt) is larger than the 
bulk etching velocity of the material {vb)- In the initial stages of the etching, 
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Figure 2: Illustration of the working principle of track etch detectors. Left: the 
production of a latent track, at the passage of a particle losing energy above the 
detection threshold. Middle: the etched cones obtained by chemical etching. 
Right: The holes produced after a longer etching. 

two cones are formed on both sides of the detector. The geometry of those cones 
depends on the REL of the incident particle; the measurement of the base area 
or of the length of the cones allows, trough a suited calibration, to determine it. 
If the etching is prolonged, the two cones form a hole in the material; the REL 
information is lost, but the detector becomes more appropriate for fast scanning 
operations. 

NTDs may be calibrated using beams of relativistic ions. A typical set-up 
consists of few sheets of detectors upstream of some target used for the beam 
fragmentation, and a set of downstream sheets of NTDs that would record the 
tracks of the non-fragmented beam ions as well as of the lighter ions produced 
in the target. Fig. 3 shows the distribution of the averaged areas of the etch 
cones in CR39 (averages are made on only two detector faces) produced by 
158 AGeV In^^ and its fragments. The exposure was made at the CERN SPS. 
By computing the REL values coresponding to different relativistic ions, the 
calibration curves for NTDs are obtained Such curves, for the CR39 (the 
squares) and for the MakrofoP (the circles) used in SLIM are presented in Fig. 
4. The variable on the ordinate refers to the so called "reduced etching rate", 
p = vt/vb- The thresholds of the two detectors corespond to p = 1. 

The strategy to search for exotic candidates in SLIM is the following: firstly, 
we perform a "strong" chemical etching of the upper CR39 detector in each 
module. In those conditions, large tracks (holes) are obtained, allowing an easy 
fast optical scan of the entire detector surface. If tracks are found, the other 
two CR39 sheets are "softly" etched (in order to obtain measurable etched 
cones along the possible track) and scanned in the areas predicted by the track 
in the first sheet. In the presence of tracks, the REL values and the orientation 
of tracks are measured. In order to accept a candidate, double coincidences 
(between the detectors above and bellow the Al absorber) are requested. If this 
is the case, the Makrofol foils are etched and scanned too. During the tests 

■^This is the first calibration of Makrofol based on the base-cone areas [2]. 
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Figure 3: Base- area distribution of the etch cones produced in the CR39 NTD 
(averaged over two detector faces) produced by Indium ions at 158 AGeV and 
their fragments. 



already performed on part of the SLIM modules, no such coincidences were 
found. 



3 Strange Quark Matter in the cosmic rays: stran- 
gelets and nuclearites 

SQM could be the ground state of quantum chromodinamics |10j . It is assumed 
that SQM is made of u, d and s quarks in nearly equal proportions. As the 
chemical potential of the s quarks in SQM is slightly larger than for u and d 
quarks, SQM is always positively charged, so electrons could neutralize it. For 
small SQM nuggets (M < 10^ GeV) the electrons would form an electronic 
cloud around the quark core; for larger masses some electrons or, for a quark 
bag radius R > lA; M > 8.4 x lO^GeV [TTl, all the electrons would be in 
equilibrium inside the SQM [El [SI 

SQM is expected to have a density slightly larger than ordinary nuclear 
matter ^1 the relation between the mass M of SQM lumps and their 
baryonic number A would be (GeV) < 0.93^. 

It was hypothesized that nuggets of SQM, with masses from those of heavy 
nuclei (in this mass region we are going to call them strangelets) to much higher 
values {nuclearites), produced in the Early Universe or in violent astrophysical 
processes, could be present in the cosmic radiation 

An upper limit for the flux of nuclearites may be obtained assuming that 
they represent the main contribution to the local Dark Matter (DM) density, 
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Figure 4: The callibration curves for the CR39 (the squares) and Makrofol (the 
circles) NTDs used in SLIM. 
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where v and M are the nuclearite average velocity and mass, respectively. 

Calculations describing the production (through binary strange stars tidal 
disruption) and the galactic propagation of cosmic ray nuclearites were recently 
published |15| . The results could be valid as orders of magnitude for the entire 
mass range of interest. 



3.1 Strangelets 

SQM should be stable for all masses larger than about 300 GeV TT. Strangelets 
with masses up to at least the multi-TeV region could be ionized and could be 
accelerated to relativistic velocities by the same astrophysical mechanisms of 
normal nuclei of the primary CR. 

They would interact with detectors (in particular NTDs) in ways similar to 
heavy ions, but with different Z/A. In ref. 12 SQM is described in analogy 
with the liquid-drop model of normal nuclei; the obtained charge versus mass 
relation is shown in Fig. 5A by the solid line, labeled "(1)". Other authors 
found different relations: Z ~ 0.1^1 for A < 700 and Z ~ 8^^/^ for larger 
masses this charge to mass relation is shown in Fig. 5 as the dashed line, 
labeled "(2)". In fT5 it was assumed that quarks with different color and flavor 
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Figure 5: A: Strangelet electric charge versus mass for different hypotheses 
discussed in refs. |E1 See text for details. B: Expected fluxes for 

strangelets in the CR near the Earth. The solid line corresponds to the assump- 
tion that their abundances follow the same rule as heavy CR nuclei The 
dashed and dot-dashed lines are from Ref . j and refer to "normal" and CFL 
strangelets. 

quantum numbers form Cooper pairs inside the SQM, the so-called color-flavor 
locked (CFL) phase, increasing the stability of the strangelets. In this case, the 
charge relation would he Z ~ O.SA^^^, shown as the dash-dotted line in Fig. 5A 
labeled "(3)". 

Several CR experiments reported possible candidate events that would sug- 
gest anomalously low charge to mass (Z/A) ratios, which could correspond to 
those expected for SQM J2] • Such candidates are reviewed in |17l IT^ . As 
strangelets could have the same origin as CR heavy nuclei, their abundances in 
the cosmic radiation could follow the same mass dependence, $ oc M~^-^, [TT]. 
The existing candidates do not contradict such an hypothesis. The solid line 
in Fig. 6 is the expected flux versus strangelet mass, assuming that the above 
assumptions are correct. 

Different nuclearite flux estimates were recently published 1151. They are 
based on the hypothesis that large nuclearites (with masses 10~^ — 10~^ so- 
lar masses) are produced in binary strange stars systems, before their gravita- 
tional collapse. The propagation inside the galaxy considers also the escape, 
spallation (through which smaller nuclearites are produced) and re-acceleration 
mechanisms. Nuclearite decays are not considered, as SQM is supposed to be ab- 
solutely stable. The predicted strangelet fluxes around the Earth are presented 
in Fig. 5B for "normal" and CFL strangelets as the dashed and the dot-dashed 
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lines, respectively. The differences originate from the different charge-to-mass 
ratios. 

The CR39 used in SLIM is sensitive (in the conditions of the "strong" etch- 
ing) to particles with REL > 200 MeV g~^ cm^ ^. This implies a minimum 
A (at the level of the detector) between 200 and 600, depending on the chosen 
A/Z model. 

If strangelets would interact with the Earth's atmosphere in the same way 
as CR nuclei, they would not reach experiments at mountain altitude. Different 
theoretical scenarios, both based on the SQM stability, were introduced in order 
to allow their deep penetration in the atmosphere; none of those mechanisms 
would allow them anyway to reach sea level. 

Mass and size decrease of strangelets during propagation. In ^7] it 
was assumed that strangelets could penetrate the atmosphere if their size and 
mass are reduced through successive interactions with the atomic air nuclei. 
This scenario is based on the spectator-participant picture. Two interaction 
models are considered: quark-quark ("standard"), and collective ("tube-like"). 
At each interaction the strangelet mass is reduced by about the mass of a Nitro- 
gen nucleus (in the "standard" model), or by more (in the "tube- like model"), 
while the spectator quarks form a lighter strangelet that continues its flight with 
essentially the same velocity. Once a critical mass is reached {A ~ 300 - 400) 
neutrons would start to evaporate from strangelets; for A < 230 the SQM would 
become unstable and decay into normal matter. In ref. |19| an estimate was 
made of the sensitivity of a NTD experiment at Chacaltaya: the mass number of 
a nuclearite penetrating the atmosphere down to that altitude would be reduced 
by a factor ~ 1/7. 

In this scenario, the minimum strangelet A value at the top of the atmosphere 
should be between 1400 and 4200, in order to be detected in SLIM. Assuming 
different strangelet structure and flux hypothesys, the expected fluxes in the 
Earth's vicinity would range between some 10~^^ and 10~^^ cm~^s~^sr^^. More 
detailed calculations are given find in Ref. Due to its expected sensitivity, 
SLIM may discern between different combinations of flux - structure hypotheses. 

Accretion of neutrons and protons during propagation. A com- 
pletely different propagation scenario was proposed in . The authors assume 
that strangelets would pick-up nuclear matter during interactions with air nu- 
clei. After each interaction, the strangelet mass would increase by about the 
atomic mass of Nitrogen, with a corresponding reduction of velocity. As the 
mass grows larger, the loss in velocity becomes smaller. They estimate that a 
strangelet of an initial A ~ 64 and an electric charge of about +2 could arrive 
at about 3600 m a.s.l. with A 340 (3600 m is the altitude of a proposed NTD 
experiment in Sandakphu, India |18p. This mechanism would also imply an 
increase of the electric charge of the strangelet, thus an increase of the Coulomb 
barrier; this may be its main difficulty. The stability of low mass strangelets 
is another questionable aspect of such a model; the expected fiuxes would be 

•^The threshold of a NTD depends on the etching conditions. A relatively high threshold for 
CR39 was chosen in order to reduce the background due to recoil tracks, neutron interactions 
and the ambient radon radioactivity. 
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larger than in the fragmentation scenario, but, for this reason, they are hard 
to estimate. Such large fluxes seem to be in disagreement with the preliminary 
SLIM results, presented in the last Section of this paper. 



3.2 Nuclearites 

In JI] was postulated that elastic collisions with atoms and molecules of the 
traversed medium are the only relevant energy loss mechanism of non-relativistic 
nuclearites with large masses, 

where p is the density of the traversed medium, v is the nuclearite velocity and 
a is its cross section: 

7r(3Af/47rpAr)2/3 for M > 8.4 x lO^^ GeV (corresponding to i?Ar ~ 1 A) 
TT X lO^^^cm^ for lower mass nuclearites 

(3) 

with pn = 3.6 X 10"'^^ g cm The cross section for nuclearites with M < 8.4 x 
10^'' GeV is determined by their electronic cloud. For nuclearites gravitationally 
trapped in our Galaxy, the average velocity would be /3 = v/c ~ lO^'^. 

Very large nuclearites (M > 3 x 10^^ GeV) could cross the Earth; as the 
SLIM sensitivity is above the DM flux limit (Eq. ^) we consider here only 
down-going nuclearites. 

A nuclearite of mass M entering the atmosphere with an initial velocity 
vq <C c, after crossing a depth L will be slowed down to 

v{L) = voe-^ ^o'' P(=')dx 

where p{x) is the air density. In the following we consider the parametrization 
of the standard atmosphere from PH). 

p{h) = ae~~ = ae ~ , (5) 

where h is the altitude, a = 1.2 x 10~'^g cm~'^ and b ~ 8.57 x 10^ cm; H is the 
total height of the atmosphere ( ~ 50 km). The integral in Eq. ^may be solved 
analytically. 

Fig. 6A shows the velocity with which nuclearites of different masses reach 
heights corresponding to typical balloon experiments (~ 40 km), to SLIM, (5.29 
km) and sea level. A computation valid for MACRO |2j (at a depth of 3400 
mwe) is also included. The velocity thresholds for detection in CR39 and in 
Makrofol are shown as the dashed curves. The decrease of the velocity thresholds 
for nuclearite masses larger than 8.4 x 10^^ GeV is due to the change in the 
nuclearite cross section, according to Eq. An experiment at the Chacaltaya 
altitude lowers the minimum detectable nuclearite mass by a factor of about 2 
with respect to an experiment performed at sea level. If the mass abundance of 
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Figure 6: A: Solid lines: arrival velocities of IMNs at different depths versus 
nuclearite mass, assuming an initial velocity outside the atmosphere of /3 = 
10~^. The nuclearites are supposed to come from above, close to the vertical 
direction. The dashed lines show the detection thresholds in CR39 (in the SLIM 
etching and Makrofol. B: Nuclearite detection conditions in CR39 (solid curves) 
and Makrofol (dashed curves), for experiments located at different altitudes. 



nuclearites decreases strongly with increasing mass this could yield an important 
increase in sensitivity. 

More general nuclearite detection conditions in CR39 and Makrofol (ex- 
pressed as the minimum velocity at the top of the atmosphere versus the nu- 
clearite mass) for different experimental locations are shown in Fig. 6B. In 
this case, the constraint is that nuclearites have the minimum velocity at the 
detector level in order to produce a track**. 

Searches for nuclearites (mostly IMNs) were performed by different experi- 
ments I23E21- The best flux upper limit was set by the MACRO experiment: 
for nuclearites with (3 ~ IQ-^ and lO" GeV < M < 10^2 GeV, the 90% C.L. 
upper limit is at the level of 2 x 10~^^ cm~2sr^^s^* , as a byproduct of the 
search for GUT magnetic monopoles |23l I24| . 

4 Preliminary results and conclusions 

Till now, we analyzed about 214 m^ of the SLIM detector, with an averaged 
exposure time at Chacaltaya of 3.5 years. As no candidate was found, the 
present 90% CL upper limit for a flux of downgoing strangelets and nuclearites 

*The CR39 threshold in the MACRO case was lower than for SLIM, due to the very low 
background at Gran Sasso Lab. 
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with M > 3 X 10^^ GeV, valid also for relativistic monopoles, is 3 x 10~^^ 
cm~^s~^sr~^. This limit disfavors the hypothesis of the accretion of matter by 
strangelets going down in the atmosphere, and can constrain some production 
or propagation models. 

We intend to complete the analysis by the end of 2006. Even if no magnetic 
monopole or SQM candidate will be found, SLIM will yield significant limits 
in mass regions not yet explored by other experiments, and will impose strong 
constraints on different scenarios describing the production and propagation of 
strangelets. 

We aknowledge many useful discussions with other members of the SLIM Col- 
laboration. Special thanks are due to the technical staff of the NTD lab. of INFN 
Bologna. 
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